Surface critical behavior of semi-infinite systems with cubic anisotropy at the 

ordinary transition. 
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^ , The critical behavior at the ordinary transition in semi-infinite n-component anisotropic cubic 

^ ■ models is investigated by applying the field theoretic approach in d = 3 dimensions up to the two- 

' loop approximation. Numerical estimates of the resulting two-loop series expansions for the critical 

exponents of the ordinary transition are computed by means of Pade resummation techniques. For 
n < He the system belongs to the universality class of the isotropic n-component model, while for 

Qn > He the cubic fixed point becomes stable, where ric < 3 is the marginal spin dimensionality of 
I the cubic model. The obtained results indicate that the surface critical behavior of the semi-infinite 

. systems with cubic anisotropy is characterized by a new set of surface critical exponents for n > n^. 

PACS numbers: PACS number(s): 64.60.Fr, 05.70.Jk, 68.35.Rh, 75.40.Cx 

O ' I. INTRODUCTION 

, The investigation of the critical behavior of real systems is an important task of the condensed matter theory. The 
' critical behavior in such systems as polymers, easy-axis ferromagnets, superconductors, as well as superfluid ^He, 
I [ Heisenberg ferromagnets, and quark-gluon plasma is described by isotropic 0{n) model vifith n = 0,1,2,3 and 4, 
' O ■ respectively, and has been analyzed in the framework of different theoretical and numerical approaches. 
^ ' Investigation of the critical behavior of real cubic crystals is one of the topics of extensive theoretical work during 
the last three decades. In crystals, due to their crystalline structure, some kind of anisotropy is always present. One 
of the simplest examples is the cubic anisotropy. A typical model of the critical behavior of such systems is the model 
with a cubic term ff added to the usual 0{n) symmetric ^{J27=i I P)^ term^-^'^. This n-component 

cubic model is a particular case of an mn-component model^, with the cubic anisotropy at m ~ 1. The model exhibits 
several types of continuous and first-order phase transitions depending on the number of spin components n, space 
dimensionality d, and the sign of the cubic coupling constant vq. The cubic models are widely applied to the study of 
magnetic and structural phase transitions. In the limiting case of rt ^ 0, it describes the critical behavior of random 
Ising-like systems^ . The case m — and n = formally describes the critical behavior of long flexible polymer chains 
_ in good solvents as model of self-avoiding walks (SAW) on a regular lattice, with short range correlated quenched 
. disorder. As has been shown by Harris^ and Kim^, the short range correlated (or random uncorrelated pointlike) 
disorder is irrelevant for such model. The case m — 1 and n ^ oo corresponds to the Ising model with equilibrium 
magnetic impurities®. 

Depending on the sign of the cubic coupling constant vq, the two types of order are possible: along the diagonals 
the type of a hypercube in n dimensions for uq > or along the easy axes of the type [1,0,...0] for vq < 0. 

' In the latter case the system can undergo the first order phase transition, as was confirmed in experiments^. In the 
' present work we are concerned with the case vq > 0. 
O ■ The presence of a surface leads to the appearance of additional complications. The source of these problems is 
^ ' connected with both the loss of translational invariance and the presence of a boundaries. General reviews on surface 
^ : critical phenomena are given in Refs.^^'^^^^^. The simplest model of critical phenomena in systems with a single planar 
surface is the semi-infinite model^°. As it is known^°'^^'^^, the phase diagram of such model is richer than that of 
^ its bulk correspondant. In general case of pure semi-infinite model with continuous 0{n) symmetries, there are a 
surface- and bulk-disordered phases (SD and BD, respectively), as well as either the surface-ordered, bulk-disordered 
phase (SO and BD, respectively) and a surface-ordered, bulk-ordered phase (SO and BO). The surface phase can 
actually occur, if d > 2 and n = 1 or d > 3 and n > 2. The boundaries between the phases are the lines of surface, 
ordinary, and extraordinary transitions which meet at a multicritical point (raQ,co) = (toqj,, c*p), representing the 
special transition and called the special point. Each of the above mentioned transitions is characterized by own fixed 
point. The constant cq is related to the surface enhancement, which measures the enhancement of the interactions at 
the surface. The coupling mg is defined in Eq.(2.1) below. The nonexistence of this ordered surface phase assume the 
occurrence of extraordinary transitions. In the case n = 2 and c? = 3 a surface transitions of the Kosterlitz-Thouless 
type is present. We do not consider this type of transition and extraordinary ones, because they have a different 
nature than the special and ordinary transitions. 

In general, there are different surface universality classes, defining the critical behavior in the vicinity of the 
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system boundaries, at temperatures close to the bulk critical point (t = (T — Tc)/Tc 0). Each bulk universality 
class divides into several distinct surface universality classes. Three surface universality classes, called respectively 
ordinary (cq — ^ oo), special (cq — c*^) and extraordinary (cq — > — oo), are known^^'^^'^^. 

In order to investigate the critical behavior of real cubic crystals we must take into account that two types of 
anisotropy can be present for such systems. The first one is a bulk anisotropy, which can be included into the 
consideration with a help of above mentioned cubic term. The other one is a surface anisotropy which arise as a 
consequence of presence the bulk cubic anisotropy (see Appendix A for details). In the present paper we are interested 
in investigation of the critical behavior only at the ordinary transition, where the surface orders simultaneously with 
the bulk. In this case, as was founded by Diehl and Eisenriegler^''''^^, the surface anisotropy is irrelevant. 

Theory of critical behavior of individual surface universality classes is very well developed for pure isotropic 
systems^^'^*'^'''^^'^^'^^, the systems with quenched surface-enhancement disorder^^'^^'^'^ and the systems with a random 
quenched bulk disorder for both the ordinary and the special surface transitions^^'^'"^. General irrelevance-relevance 
criteria of the Harris type for the systems with quenched short-range correlated surface-bond disorder were predicted 
in^^ and confirmed by Monte-Carlo calculations^^'^^. Moreover, it was established that the surface critical behavior 
of semi-infinite systems with quenched bulk disorder is characterized by the new set of surface critical exponents in 
comparison with the case of pure systems^'''^^. 

The remainder of this paper is organized as follows. Section II contains the description of model and further useful 
background. In Section III the renormalization group approach is described. Section IV contains the calculations of the 
surface renormalization factor Zg^ and surface critical exponent r]Q^ by applying the field theoretic approach directly 
in = 3 dimensions, up to the two-loop order. The numerical estimates of the resulting two-loop series expansions 
for the critical exponents of the ordinary transition are presented in Section V. The calculations are performed by 
means of the Pade resummation techniques for the cases n = 3, 4, 8 and for the case of n — > oo, which corresponds 
to the Ising model with equilibrium magnetic impurities. Section VI contains concluding remarks. The Appendix 
A contains the Landau-Ginzburg- Wilson functional for a semi-infinite systems of spins with cubic anisotropy. The 
Appendix B contains the standard surface scaling relations for the case d = Z. 



II. THE MODEL 



The efii'ective Landau-Ginzburg- Wilson Hamiltonian of the n- vector model with cubic anisotropy in the semi-infinite 
space is given by (see Appendix A) 

oo 

H{$) = jdzj d'-\[\ I V^p +\ml \$\^ 



+ -^v<yY.\^it + ^MY.\c^i?f] (2.1) 

i=l ■ i=l 

where (j){x) = {^i(a;)} is an n-vector field with the components (j>i{x), i = l,...,n. Here mg is the "bare mass", 
representing a linear measure of the temperature difference from the critical point value. The parameters uq and vo 
are the usual "bare" coupling constants uq > and > 0. It should be mentioned that the d-dimensional spatial 
integration is extended over a half-space IR^^ = {x=(r, z) € M'^, with r e 1R'^~^ and z > 0}, bounded by a planar 
free surface at ^; = 0. The fields (/>i(r, z) satisfy the Dirichlet boundary condition <pi{r, 2) = at = in the case of 
ordinary transition, and the Neumann boundary condition dn(f>i(r, z) = at z = in the case of special transition^**' 
The model defined in (2.1) is translationally invariant in directions parallel to the external surface, z = 0. Thus, we 
shall use mixed representation, i.e. Fourier representation in d — 1 dimensions and real-space representation in z 
direction. Therefore, as was mentioned above*''''*® and shown in Appendix A, no specific surface term will appear in 
the case of ordinary transition, when the Dirichlet boundary conditions on the surface are assumed. 

The added c;ubic; term breaks the 0(n) invariance of the model, leaving a discrete cubic symmetry. The model 
(2.1) has four fixed points: the trivial Gaussian, the Ising one in which the n components are decoupled, the isotropic 
(O(n)-symmetric), and the cubic fixed points. The Gaussian and Ising fixed points are never stable for any number 
of components n. For isotropic systems, the 0(n)-symmetric fixed point is stable for n < Uc, whereas for n > ric it 
becomes unstable. Here Uc is the marginal spin dimensionality of the cubic model, at which the isotropic and cubic 
fixed points change stability, i.e. for n > ric, the cubic fixed point becomes stable. The 0(n)-symmetric fixed point is 
tricritical. At n = Uc, the two fixed points should coincide, and logarithmic corrections to the 0(n)-symmctric critical 
exponents are present. The calculation of the critical marginal spin dimensionality is the crucial point in studing 
the critical behavior in three-dimensional cubic crystals. Diflferent results for ric have been published in a series of 
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works in which different methods have been used. In the framework of the field-theoretical RG analysis the one-loop 
and three- loop approximations at e = 1 lead to the conclusion that should lie between 3 and 4^'^''^^, and the cubic 
ferromagnets are described by the Heisenberg model. On the other hand, by using the field theoretic approach directly 
in d = 3 dimensions up to the three-loop approximation, it has been found that rig = 2.9^^'^^. Similar conclusions were 
obtained in^'^, where it was found that Uc = 2.3. The calculations performed by Newman and Riedel^* with the help 
of the scaling-field method, developed by Goldner and Ridel for Wilson's exac;t momentum-space RG equations, have 
given for d = 3 the value ric = 3.4. Field-theoretical analysis, based on the four-loop series in three dimensions'^'""' 
and results of the five- loop e = 4 — expansion^°'^^'^^ suggest that ric < 3. Recently, a very precise six-loop result 
for the marginal spin dimensionality of the cubic model, ric = 2.89(4), was obtained in the framework of the 3D 
field-theoretic approach**'. Thus, it was finally established that the critical behavior of the cubic ferromagnets is not 
described by the isotropic Heisenberg Hamiltonian, but by the cubic model, at the cubic fixed point. However, it was 
found that the difference between the values of the bulk critical exponents at the cubic and the isotropic fixed points 
is very small, i.e. it is hard to determine this difference experimentally. Nevertheless, the recently obtained results 
stimulated us to perform the analysis of the surface critical behavior of semi-infinite n-component anisotropic cubic 
model, and to determine corresponding surface critical exponents. 



III. RENORMALIZATION 



The fundamental two-point correlation function of the free static theory corresponding to (2.1) is defined by the 
Dirichlet propagator: 

{(Pi{r, z)ipj{0, z'))o = Gnir; z, z')6ij. (3.1) 
In the mixed pz representation the Dirichlet propagator is 

Gd{p; z, z') = L---o\^-^'\ - e-«o(z-K^')l ^ (3.2) 
2ko L J 

where the standard notationis used and, kq = \/p'^ + rn^. The propagator vanishes identically when at least one of its 
z coordinates is zero, because we have assumed the Dirichlet boundary conditions. Consequently, all the correlation 
functions involving at least one field at the surface vanish. This property holds for both the free and the renormalized 

theories^^. 

In fact the critical surface singularities at the ordinary transition can be extracted by studying the nontrivial in 
this case correlation function involving the {inner) normal derivatives of the fields at the boundary, 0„^(r)^^'^^'^^. 
Actually, in order to obtain the characteristic exponent rj°^''^ of surface correlations, it is sufficient to consider a 
correlation function with two normal derivatives of boundary fields, i.e. 

^2(rt = ^ij>,z)l^, A ^(-p,.')L.=o) , (3-3) 

where the fields </?(p, z) are the Fourier transforms of the fields iyj(r, z) in the {d — 1) dimensional parallel to surface 
space. Q2{p) is a parallel Fourier transform of the corresponding two-point function Q2{t) in direct space. At 

the critical point Q2{p) behaves as . It reproduces the leading critical behavior of a two-point function 

G2{p) = {fip, z)(p{—p, z )) in the vicinity of the boundary plane. The surface critical exponent 77^'"'^ is provided by 
the scaling dimension of the boundary operator 9„(y?(r). 

The surface correlation function exponent 77^'"'^ in the semi-infinite systems with cubic anisotropy differs from its 
corresponding value for the isotropic semi-infinite system. The remaining surface critical exponents of the ordinary 
transition can be determined through the surface scaling laws*^. (see Appendix B). 

In the present formulation of the problem, the renormalization process for the cubic anisotropic system is essentially 
the same as in the isotropic casc*^'^°. Explicitly, the renormalized bulk field and its normal derivative at the surface 
should be reparametrized by different uv-finite renormalization factors Z^(u,v) and Zg^p{u,v) 

(Pr{x)=Z~^ <p{x) and {d„ip{r)) j^=Z~^ d„tp{r), 
and renormalized correlation functions involving N bulk fields and M normal derivatives are 

g\i"'''''\p;m,u,v) = Z-^Z-^g^^'^\p;mo,uo,vo) (3.4) 
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for {N,M) 7^ (0,2). In order to remove the ultraviolet (uv) singularities of the correlation function Q^^'^^ with 
two surface operators {N, M) = (0, 2) in the vicinity of the surface, an additional, additive renormalization (zero- 
momentum subtraction) is required, so that 



7(0,2) 



(p)-e(°'')(p=o) 



(3.5) 



The typical bulk uv singularities, which are present in the correlation function are subtracted via the standard 

mass renormalization of the massive infinite-volume theory. It also relates to coupling constants, for which standard 
vertex renormalization of coupling constants takes place. 

The siuface renormalization factor Zg^{u, v) can be conveniently obtained from the consideration of the boundary 
two-point function 5^"'^^ 



TO d 



P^o p dp 



(3.6) 



A standard RG argument involving an inhomogeneous Callan-Symanzik equation yields the anomalous dimension of 
the operator 9„<^(r) 



d 

= A(n,.) ^^"^f^("'"^ + /?.(n,.) ^^"^^^^"'") 



(3.7) 



du 



dv 



FP 



"FP" indicates here that the above value should be calculated at the infrared-stable cubic fixed point of the underlying 
bulk theory, {u,v) = {u*,v*). The surface critical exponent r/^'"'* at the ordinary transition is then given by 



7?f<^ = 2 + r?a^. 



(3.8) 



IV. PERTURBATION THEORY UP TO TWO-LOOP APPROXIMATION 

After performing the mass and additive renormalization of the correlation function tj*^"'^^ (p) and carrying out the 
integration of Feynman integrals by analogy with Ref.^"'^'*, we obtain for renormalization factor 



j{0) 

Z9Uuo,vo) = 1 + ^ + H^'C, (4.1) 



where the constant C follows from the two-loop (melon-like diagrams) contribution to the correlation function and 
has the value, 

107 7 4 

C~ — - -In- - 0.094299 ~ -0.105063. (4.2) 
162 3 3 ^ ^ 

The coefficients 4°'' and are the weighting factors belonging to one- and two-loop (melon-like) diagrams in the 
Feynman diagrammatic expansion of the correlation function g^'^''^\p), and equal 

- ^, with t^^^ = 'l±luo + vo, (4.3) 
with 4°^ = ul + vl + 2voUo ■ (4.4) 

The factors a^iid 4°^ follow from the standard symmetry properties of Hamiltonian (2.1). Here the renor- 
malization factor Zg^ is expressed as a second-order series expansion in powers of bare dimensionless parameters 
uo = uo/{8Trm) and vq = vo/{8i:m). As it is usual in swperxenormalizable theories, the renormalization factor 
expressed in terms of unrenormalized coupling constants, is finite. 
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As a next step, the vertex renormalizations should be carried out. To the present accuracy, they are 

n + 8 



Wo = u\l-\ —u + vy (4.5) 

= v{l+'^v + 2v}j . (4.6) 



As known, the vertex renormalization at li = 3 is a finite reparametrization. All relevant singularities have been re- 
moved already after the mass renormalization and taking into account the special bubble-graph combinations emerging 
in the theory with Dirichlet propagators. Thus we obtain a modified series expansion up to two-loop approximation 



ry -, n+2 _ V n+2 / n+8\ 



+ {C+-j v' + 2^C+—j uv. (4.7) 

Combining the renormalization factor Zgip{u, v) together with the one- loop pieces of the beta fmictions : Pu{u, v) = 
— w^l — u — ^ and I3y{u, v) = —v(^ — — 2v}j and inserting them into Eq.(3.7), we obtain the desired series 
expansion for r]Q^, 

n-\-2 V 



2(n+S) 6 

ti!^C(n)w^ - ^C(l).^ - 
(n+8)2 ^ ' Q ^ ' n+8 



where C{n) is a function of the order-parameter components number n, and is defined as 

C{n)=C+-^, (4.9) 

whereas the renormalized coupling constants u and v, normalzed in a standard fashion are w=^^^w and f =f w- 

Eq.(4.8) supplies our result for the critical exponent of surface correlation function for the model with the effective 
Hamiltonian of the Landau-Ginzburg- Wilson type with cubic anisotropy in the semi-infinite space (2.1) with general 
number n of order parameter components. 

The knowledge of r]\\ gives possibility to calculate the other surface critical exponents through the scaling relations. 
For convenience, from now on we omit the superscript ord for the surface critical exponents. 

The critical exponents should be calculated for different n (n = 3, 4, 8, and n oc) at the standard infrared-stable 
cubic fixed (FP) points of the underlying bulk theory, as is usually accepted in the massive field theory. As was 
mentioned above, in the cases n < ric the cubic ferromagnets are described by the Heisenberg isotropic Hamiltonian 
at the 0(n)-symmetric fixed point. 

In the case of the replica limit n ^ we obtain from (4.8) the series expansion of ryJi" for semi-infinite random 

Ising-like systems. This case was investigated in detail of one of us previously^^. 



V. NUMERICAL RESULTS 



In order to obtain the full set of surface critical exponents for the ordinary transition in systems with cubic 
anisotropy, we substitute the expansion (4.8) for t^h into the standart scaling-law expressions for the surface exponents 
(see Appendix 1). 

For each of the above mentioned surface critical exponents of the ordinary transition we obtain for d = 3 a double 
series expansion in powers of u and w, truncated at the second order. As it is known^*'^^'®*^'^^ , power series expansions 
of this kind are generally divergent due to a nearly factorial growth of expansion coefficients at large orders of 
perturbation theory. In order to perform the analysis of these perturbative series expansions and obtain accurate 
estimates of the surface critical exponents, a powerful rcsummation procedure must be used. One of the simplest 
ways is to perform the double Pade-analysis"'^. This should work well when the series behaves in lowest orders "in a 
convergent fashion" . 

The results of our calculations of the surface critical exponents of the ordinary transition for various values of 
n = 3,4, 8, (X) at the corresponding cubic fixed points are presented in Tables 1-5. Unfortunately, the second-order 
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(p=2) analysis of pcrturbativc scrics^"^ gives the cubic fixed point with coordinates uq = 1.5347 and vq = —0.0674 
at ri = 3 for the 3D model. The analysis of the eigenvalues of the stability matrix shows that in the frames of the 
two-loop approximation the cubic fixed point at n = 3 is unstable and the 0(n)-symctric fixed point is stable. But, 
the estimates of the marginal spin dimensionality of the cubic model Uc in the frames of three- loop^^'^^, four-loop^^'^°, 
five-loop e = 4 — d expansion*^'^^'^^ and six- loop study show that the cubic ferromagnets are not described by the 
Heisenberg isotropic model, but by the cubic model at the stable cubic fixed point. Higher precision six-loop field- 
theoretical analysis*^ give the value of the marginal spin dimensionality of the cubic model equal to ric = 2.89(4). 
In accordance with this we use the cubic fixed point of the higher p = 3 order of perturbative series for obtaining 
the set of surface critical exponents at n = 3. For estimation of the reliability of the obtained results we performed 
calculations at the cubic fixed point of the p = 6 order in Table 2. We obtained that difference in these two cases 
are approximatively: 0.5% for r]\\, 0.4% for r]^ 0.8% for Ai, 0.1% for 6% for 7n, 0.2% for 71, 0.2% for 5i, and 
0.8% for 5ii. The obtained results indicate that the difference in the ways of the (3 fimctions resummation have no 
essential influence on the values of the surface critical exponents and that the results obtained in the frames of the 
two-loop approximation are stable and reliable. The surface critical exponents of the ordinary transition for n = 4, 8 
and n — »■ 00 were calculated at the standard infrared-stable cubic fixed (FP) points of the underlying bulk theory, as 
it is usually accepted in the massive field theory. 

The quantities O1/O2 and OujO^i represent the ratios of magnitudes of first-order and second-order perturbative 
corrections appearing in direct and inverse series expansions. The larger (absolute) value of these ratios indicate the 
better apparent convergence of truncated series. 

The values [p/q] (where p, g = 0, 1) arc simply Fade approximants which represent the partial sums of the direct and 
inverse series expansions up to the first and the second order. The nearly diagonal two- variable rational approximants 
of the types [11/1] and [1/11] give at u = or u = the usual [1/1] Fade approximant*^. As it is easy to see 
from Tables 1-5, the values of [11/1] and [1/11] Fade aproximants do not differ significantly between themselves. Wc 
consider these values as the best we could achieve from the available knowledge about the series expansions in the 
two-loop approximation scheme. Thus, our final results are presented in the last columns of the Tables 1-5. Their 
deviations from the other second-order estimates might serve as a rough measure of the achieved numerical accuracy. 
As it is easy to see, the obtained results indicate about good stability of the results calculated in the frames of the 
two- loop approximation scheme. 

The results for surface critical exponents of semi-infinite model with cubic anisotropy, calculated at the cubic 
fixed point are different from the results for surface critical exponents of standart semi-infinite n-component model 

(see".17,20,57)^ 

If n < ric, the cubic fixed point is unstable and the cubic term in the Hamiltonian (2.1) becomes irrelevant. In this 
case the isotropic fixed point is stable and the system is described by the simple 0(n)-symmetric model in 3D. The 
corresponding surface critical exponents can be calculated from the series, presented in^". 

As was indicated previously, in the limit n ^ 0, the cubic model (2.1) with mq < and iiq > describes the 
semi-infinite Ising-like systems with random bulk disorder. The investigation of the ordinary transition for such kind 
of systems was presented in^**. 

VI. CONCLUDING REMARKS 

We have studied ordinary transition for a semi-infinite systems with cubic anisotropy by applying the field theoretic 
approach directly in d = 3 dimensions, up to the two-loop approximation. We have performed a double Fade analysis 
of the resulting perturbation series for the surface critical exponents of the ordinary transition for various n = 3, 4, 8, 00, 
in order to find the best numerical estimates. Wc find that at n > nc the surface critical exponents of the ordinary 
transition in semi-infinite systems with cubic anisotropy belong to the cubic universality class. 

In order to obtain more precise numerical estimates for the case of 3D dimensional cubic crystal with n = 3, a 
further theoretical investigation of the asymptotic surface critical behavior of semi-infinite cubic systems would be 
highly desirable within the framework of higher-order RG approximations. 

We suggest that the obtained results could stimulate further experimental and numerical investigations of the 
surface critical behavior of random systems and systems with cubic anisotropy. 
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Appendix A. The Landau-Ginzburg-Wilson functional for the semi-infinite systems with cubic anisotropy 

We outline the main steps of deriving the Landau-Ginzburg-Wilson functional for the case of Heisenberg ferromagnet 

and discuss briefly the surface term leading to the surface boundary conditions. We perform the calculation for a 
simple cubic structure and start from the mean-field approximation for the Heisenberg exchange interaction, as the 
cubic- anisotropy term can be added afterwards. Additionally, we limit ourselves to one-component theory, as the 
generalization to the n-component version is not important to the essence of the principal argument. 

Suppose we have the system of localized spins of magnitude S and described by the Heisenberg Hamiltonian, which 
in the applied field h takes the form 



In the mean-field approximation 

Sj • Sj =< Si > •Sj+ < Sj > •Si- <Si> •<Sj >, 
and for the spin quantization axis taken as the z-axis. We can write the free energy in the form 

sinh(,'?/i,(5+l/2)) , 1 



F = -NkBTln- 



sinh(/?/ii/2) 



2^J., <5f ><5|>, 



(Al.l) 



(A1.2) 



(A1.3) 



where hi = J2j Jij < "^1 > +^ eff'ective field acting on Sf, /3 = {UbT) ^ is the inverse temperature in energy 

units, and N is the total number of spins. 

The constant term can be rewritten in the first nontrivial order the continuum-medium approximation as 



\ J2 Jij < 5f X 5| >^ 1 Jo < 5^(x) > 



-al < 5^(x) > < 5^(x) > 



x=Ri 



(A1.4) 



x=Ri 



where Jo = Jy , ao is the lattice constant, and i = Ri denotes here the lattice site position. 
The expansion of the ratio of hyperbolic function in y = (3hi can be represented as 

+ ^^(3^' + 65' + 25-l)y4] + 0(y6). 
Hence, the free energy (per site) to the same approximation after taking the continuum-medium limit reads 



(A1.5) 



F 

N 



N 



1+ (dH 
V J 



25(5 + 1 



-fcB(T-Te)<^(x)2 



+ ^Joa2|V</.(x)|2 - ^5(5+ l)/3/iJo</.(x) 



+ ^{S''{S + If - ^(35^ + 652 ^ 25 - l))(Jo.^(x)4)) 

+ 01 



S. 

^) + o{h^) + o{\V4>\% 



(A1.6) 



with ^ = — A;BTln(25 -|- 1). The continuous field ^(x) expresses the limiting value of < 5? > per volume Oq. 

Defining the mean-field critical temperature Tc = j^JciS{S + 1) and assuming that we can put T ~ in the last 
two terms we obtain the desired free energy functional in the form 



F 
N 



N 



+ ^i</,(x)r-/.</.(x)). 



(A1.7) 



where the constants are related in an obvious fashion to the coefficients TOq and mq, when we divide F by the exchange 
stiffness constant Aq, which contains non-divergent constants at the critical point. 
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Now, the boundary conditions appear when we derive the Landau- Ginzburg equation in the form of Euler equation 
for </)(x). However, in such a situation specific surface terms appear. The are two types of terms. First, is the gradient 
term ^o|V<;^(r, z = O)]"^, since in Jo = J2j Jij the spins above the surface are missing. Second, the geometrical surface 
term of the form ^cq / d'^^^r(jp{r, z — 0) may appear where the cq is the surface enhancement constant. This is 
becouse on the surface the role of the bulk cubic anisotropy is taken over by the surface anisotropy, which in the first 
nontrivial order has the form i J27=i '^o'l^ii^^^ — effect, the condition on the surface coming from the Euler 

variational scheme takes the following form 



J2 I d''-''r{A'on»y(f>i{^,z)+ci<f>i{r,z)}5ct>i{r,z)U=o = 0, (A1.8) 

i=l 

where d(j)i{r,z = 0) is the variation of (j)i{r,z) on the surface and n is the vector perpendicular to the surface. 
Therefore, we can choose the boundary conditions in either way, namely 

1) i5(/>i(r,z)|^=o = 0, i.e. 4>i{r, z)\z=o = const = So, (A1.9) 

2) (a„</)i(r, z) + -j0i(r, z))U=o = 0. (ALIO) 

In this paper we have selected the b.c. 1) with 5*0 = 0, which correspond to the Dirichlet boundary condition. In 
such a situation, no specific surface term appears in the starting functional (2.1). However, it must be said that it is 
the concrete experimental situation that determines the type of boundary conditions to be taken into the theoretical 
analysis. 



Appendix B. Scaling relations between the surface critical exponents 

The individual RG series expansions for other critical exponents can be derived through standard surface scaling 
relations^^ with d = 3: 



/3i = 




-^ + v\\), 


7ii = 






71 = 




- v±), 


Ai = 








A 


d + 2-r] 


^1 = 


Wi~ 


d-2 + r]f 


^11 = 


Ai 
" Pi ' 


_ d-v\\ 

d-2 + r]\i 



Each of these critical exponents characterizes certain properties of the cubic anisotropic system near the surface. 
The values v, t], and A = ^{d +2 — r])/2 are the standard bulk exponents. 
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TABLE I: Surface critical exponents of the ordinary transition for d = 3 up to two-loop order at the cubic fixed point (of order 
p = 3): u* = 1.348, v* = 0.074, at n = 3. 



exp 


02 


On 

0,, 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[11/1] 


[1/11] 


/ 


m 


2.27 


1.67 


2.00 


1.681 


1.725 


1.541 


1.594 


1.429 


1.428 


1.429 


v± 


2.74 


1.91 


1.00 


0.841 


0.863 


0.783 


0.805 


0.749 


0.749 


0.749 


Ai 


2.32 


3.69 


0.25 


0.409 


0.440 


0.478 


0.504 


0.530 


0.530 


0.530 


f3i 


-4.39 


-2.58 


0.75 


0.909 


0.940 


0.873 


0.858 


0.880 


0.880 


0.880 


711 


0.00 


0.00 


-0.50 


-0.50 


-0.50 


-0.433 


-0.428 






-0.400 


71 


3.38 


17.70 


0.50 


0.739 


0.814 


0.810 


0.838 


0.839 


0.837 


0.838 


<5i 


1.99 


2.53 


1.67 


1.844 


1.865 


1.933 


1.957 


2.023 


2.023 


2.023 




1.83 


2.47 


0.33 


0.475 


0.498 


0.552 


0.582 


0.647 


0.647 


0.647 



TABLE II: Surface critical exponents of the ordinary transition for d = 3 up to two-loop order at the cubic fixed point (of order 
p = 6) u* = 1.321(18), V* = 0.096(20), at n = 3. 



exp 


Oi 


On 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[11/1] 


[1/11] 


/ 


V\\ 


2.29 


1.68 


2.00 


1.684 


1.727 


1.545 


1.597 


1.436 


1.435 


1.436 


V± 


2.76 


1.92 


1.00 


0.842 


0.863 


0.785 


0.806 


0.752 


0.752 


0.752 


Ai 


2.34 


3.72 


0.25 


0.408 


0.438 


0.476 


0.501 


0.526 


0.526 


0.526 


/3i 


-4.42 


-2.60 


0.75 


0.908 


0.938 


0.872 


0.858 


0.879 


0.879 


0.879 


711 


0.00 


0.00 


-0.50 


-0.50 


-0.50 


-0.433 


-0.428 






-0.424 


71 


3.41 


17.87 


0.50 


0.737 


0.811 


0.807 


0.834 


0.835 


0.836 


0.836 


Si 


2.01 


2.55 


1.67 


1.842 


1.863 


1.930 


1.953 


2.018 


2.017 


2.018 


Sii 


1.84 


2.48 


0.33 


0.474 


0.497 


0.550 


0.579 


0.642 


0.642 


0.642 



TABLE III: Surface critical exponents of the ordinary transition for d = 3 up to two-loop order at the cubic fixed point (of 
order p = 2) m* = 1.064, v* = 0.520, at n = 4. 



exp 


Oi 

0?, 


0,, 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[11/1] 


[1/11] 


/ 


V\\ 


2.12 


1.54 


2.00 


1.647 


1.700 


1.481 


1.550 


1.319 


1.314 


1.317 


Vi- 


2.51 


1.74 


1.00 


0.824 


0.850 


0.753 


0.783 


0.705 


0.705 


0.705 


Ai 


2.07 


3.27 


0.25 


0.426 


0.464 


0.511 


0.549 


0.588 


0.589 


0.589 


/3i 


-4.95 


-2.64 


0.75 


0.926 


0.964 


0.891 


0.873 


0.898 


0.900 


0.899 


711 


0.0 


0.0 


-0.50 


-0.50 


-0.50 


-0.409 


-0.400 






0.311 


71 


2.91 


12.61 


0.50 


0.765 


0.860 


0.856 


0.900 


0.899 


0.902 


0.901 


Si 


1.83 


2.34 


1.67 


1.863 


1.889 


1.969 


2.003 


2.105 


2.103 


2.104 


Sii 


1.70 


2.31 


0.33 


0.490 


0.519 


0.583 


0.623 


0.727 


0.724 


0.726 



* Electronic address: pylyp@ph.icmp.lviv.ua 
^ A. Aharony, Phys. Rev. B 8, 4270 (1973). 

2 I.J. Ketley and D.J. Wallace, J.Phys.A 6,1667 (1973). 

3 D.J. Wallace, J. Phys. C 6, 1390 (1973). 

* E. Brezin, J. C. Le Guillou, J. Zinn- Justin, Phys. Rev. B 10, 892 (1974). 
^ G. Grinstein and A. Luther, Phys. Rev. B 13, 1329 (1976). 



10 



TABLE IV: Surface critical exponents of the ordinary transition for d = 3 up to two-loop order at the cubic fixed point (of 
order p = 2) u* = 0.525, v* = 1.146, at n = 8. 



exp 


O2 


Ol,- 

0,, 


fO/Ol 

L / J 


[1/0] 


[0/1] 

L / J 


[2/0] 

L / J 


[0/2] 


[11/1] 

L / J 


[1/11] 

I / J 


/ 








9 on 


± .U^O 


1 .uyy 


1 A7Q 


± .0^0 


1 "^07 

1 . OU 1 


1 9Q7 






2.51 


1.74 


1.00 


0.823 


0.849 


0.752 


0.781 


0.702 


0.701 


0.702 


Ai 


2.05 


3.22 


0.25 


0.427 


0.466 


0.514 


0.553 


0.592 


0.593 


0.593 




-5.45 


-2.77 


0.75 


0.927 


0.966 


0.895 


0.878 


0.903 


0.905 


0.904 


7ii 


0.0 


0.0 


-0.50 


-0.50 


-0.50 


-0.393 


-0.380 






-0.315 


71 


2.84 


11.57 


0.50 


0.766 


0.863 


0.860 


0.907 


0.902 


0.909 


0.906 


Si 


1.84 


2.34 


1.67 


1.864 


1.890 


1.971 


2.005 


2.113 


2.139 


2.126 


Sii 


1.70 


2.33 


0.33 


0.491 


0.521 


0.584 


0.624 


0.739 


0.732 


0.736 



TABLE V: Surface critical exponents of the special transition for d = 3 up to two-loop order at the cubic fixed point (of order 
p = 2) u* = 0.201, V* = 1.508, for n -» 00. 



exp 


Ol 

02 




[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[11/1] 


[1/11] 


/ 


m 


2.14 


1.56 


2.00 


1.648 


1.701 


1.484 


1.552 


1.312 


1.300 


1.306 




2.53 


1.75 


1.00 


0.824 


0.850 


0.755 


0.783 


0.706 


0.704 


0.705 


Ai 


2.07 


3.26 


0.25 


0.426 


0.463 


0.511 


0.549 


0.585 


0.587 


0.586 


/3i 


-5.34 


-2.75 


0.75 


0.926 


0.963 


0.893 


0.876 


0.901 


0.904 


0.903 


711 


0.0 


0.0 


-0.50 


-0.50 


-0.50 


-0.387 


-0.373 






-0.323 


71 


2.88 


11.93 


0.50 


0.764 


0.858 


0.856 


0.901 


0.895 


0.902 


0.899 


5i 


1.85 


2.36 


1.67 


1.862 


1.888 


1.968 


2.001 


2.108 


2.105 


2.107 


Sii 


1.71 


2.34 


0.33 


0.490 


0.519 


0.581 


0.621 


0.736 


0.728 


0.732 
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